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Negative  differential  conductivity  in  conducting  carbon  nanotubes 

Ant.  S.  Maksimenko  and  G.  Ya.  Slepyan 

Institute  for  Nuclear  Problems,  Belarus  State  University, 

Bobruiskaya  1 1 , 220050  Minsk,  Belarus 

Abstract.  A theoretical  model  and  computations  of  the  7 — V characteristics  of  long  carbon  nan- 
otubes in  a strong  axial  dc-fields  at  room  temperatures  have  been  presented.  Negative  differential 
conductivity  has  been  predicted.  It  has  been  shown  that  [d//dV|  for  metal  carbon  nanotubes  in 
the  region  of  the  negative  differential  conductivity  significantly  exceeds  corresponding  values  for 
semiconducting  ones.  The  predicted  effect  makes  possible  the  design  of  wave-generating  nanotube- 
based  diodes  for  submillimeter  and  infrared  ranges. 


Since  the  discovery  by  Iijima  [ ] of  carbon  nanotubes  (CNs).  a great  deal  of  interest  has 
been  focused  on  these  quasi-one-dimensional  monomolecular  structures  because  of  their 
unique  physical  properties  (mechanical,  electrical,  optical,  etc.)  and  the  rapid  experimental 
progress  in  the  controlled  preparation.  Processes  of  electron  transport  in  strong  external 
fields  when  nonlinear  effects  are  constitutive  are  of  a great  interest  for  potential  applications 
in  nanoelectronics  and  for  experimental  diagnostic  of  CN  themselves. 

In  this  report  we  report  a theoretical  phenomenological  analysis  of  the  /—  characteristics 
of  CNs  at  room  temperatures,  when  k#T  > £c.  A£.  In  our  case,  the  current  is  produced 
by  free  charge  carriers  — quasiparticles  which  are  tt -electrons  moving  in  the  field  of  the 
crystalline  lattice.  The  nonlinearity  of  the  I — V characteristic  appears  due  to  the  nonlinear 
properties  of  the  quasiparticle  gas.  We  investigate  the  negative  differential  conductivity 
(NDC)  d//dV  < 0 in  CNs  in  a certain  range  of  the  field  strength,  predicted  in  | ]. 

Let  us  consider  an  undoped  single-wall  achiral  CN  (m, 0)  or  (m,m)  exposed  to  a ho- 
mogeneous axial  dc-field  Ez,  Ez  = V/L,  where  V is  the  voltage  between  the  CN  ends. 
We  shall  apply  the  semi-classical  approximation  considering  the  motion  of  n -electrons  as 
a classical  motion  of  free  quasi-particles  with  dispersion  law  extracted  from  the  quantum 
theory.  With  the  account  to  the  hexagonal  crystalline  structure  of  CNs,  the  tight-binding 
approximation  gives  [ ] : 


£s(Pz)  - ±K0 
for  the  zigzag  CNs  and 
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for  the  armchair  CNs.  Here  yo  ~ 3.0  eV  is  the  overlapping  integral,  a — 3 b /2h,  b — 1.42A 
is  the  C-C  bond  length.  The  — and  + signs  correspond  to  the  valence  and  conduction  bands, 
respectively.  In  view  of  the  transverse  quantization  of  the  quasi-momentum,  its  transverse 
component  can  take  m discrete  values,  p,p  = sAp $ — nffhs/am  (,v  = 1, . . . , m)  and 
we  used  £(s  A p,/,,  pz)  = £s(pz)  in  both  equations.  As  different  from  we  assume  pz 


424 


QWR/QD.21p 


425 


continuously  varying  within  the  range  0 < pz  < 2 n/a  what  corresponds  to  the  model  of 
infinitely  long  CN  (L  — oo).  This  model  is  applicable  to  the  case  under  consideration 
because  we  are  restricted  to  temperatures  and/or  voltages  well  above  the  level  spacing  | ]. 

The  motion  of  quasi-particles  in  an  external  axial  electric  dc-field  is  described  by  the 
Boltzmann  kinetic  equation 


9 f(p)  1 

eEz^-  = —[f(p)-F(p)  1,  (3) 

3 Pz  ? 

where  e is  the  electron  charge,  F{p)  is  the  equilibrium  Fermi  distribution  function  and  r is 
the  relaxation  time.  The  relaxation  term  of  the  equation  (3)  describes  the  electron-phonon 
scattering  [ , ],  electron-electron  collisions  etc. 

Utilizing  the  method  originally  developed  in  the  theory  of  quantum  semiconductor 
superlattices  we  can  construct  an  exact  solution  of  kinetic  equation  (3)  without  assuming 
the  electric  field  to  be  weak  [ ]. 

Then  in  view  of  Eqs.  (l)-(3)  one  can  obtain  [ ] 


jz(Ez)  = 


&eY0  r2Qz 

7f^E,  + (rQT)2E^ 

r= 1 .5=1 


(4) 


Here  Q = aeEz  for  the  zigzag  CNs  and  Q = aeEz/ \/3  for  the  armchair  ones  is  the  Stark 
frequency.  This  equation  states  the  basis  for  the  evaluation  of  I — V characteristics.  Let  us 
estimate  the  restrictions  to  the  theoretical  approach  being  developed.  As  it  has  been  stated 
above,  the  quasi-particles  motion  is  described  classically  by  Boltzmann  equation  (3).  It 
imposes  the  limitation  on  the  external  electric  field  strength:  \EZ\  < yo/2eR  | :].  The 
Coulomb  electron-electron  interaction  has  been  also  left  out  of  account  in  our  analysis. 
The  role  of  this  mechanism  as  applied  to  CNs  was  discussed  in  | ]. 

The  I—V  characteristics  obtained  via  numerical  calculation  of  Eq.  (4)  are  presented 
in  the  right  panel  of  Fig.  I for  the  metal  zigzag  CNs  and  in  the  left  one  for  the  armchair 
CNs  (m  ^ 3 q)  zigzag  CNs.  The  figures  show  the  linear  dependence  of  jz  on  Ez  at  weak 
strengths  of  the  external  field;  it  corresponds  to  the  region  of  ohmic  conductivity.  As  Ez 
increases,  the  value  9 jz/dEz  growth  smaller  and  at  Ez  = Ezmx)  the  current  density  reaches 
the  maximum  value  y'7max . Further  increase  of  Ez  results  in  the  decrease  of  jz.  At  small 

m the  magnitudes  of  Ezmdx)  and  j':mx  for  the  armchair  CNs  turn  out  to  be  less  than  for 
zigzag  CNs  with  similar  radius  by  the  factor  approximately  2.  As  different  from  the  zigzag 
CNs,  the  increase  of  m leads  to  the  Ezmax>  growing  up,  and  I—V  curves  for  zigzag  and 
armchair  CNs  become  practically  identical  at  R > 20  nm.  Thus,  we  predict  the  region  with 
the  negative  differential  conductivity  djz/dEz  < 0,  in  the  I — V characteristics  of  CNs. 

The  external  field  strength  £^max)  « 3.2  • 103  V/cm  for  the  NDC  region  appears  to  be 
unexpectedly  weak  because  the  nonlinearity  in  the  structure  under  consideration  is  deter- 
mined by  the  value  of  aEz;  quantum  superlattices  with  periods  about  10~6  cm  [ |,  much 
larger  than  b , show  approximately  the  same  fields  for  the  NDC  manifestation.  It  means 
that  the  nonlinearity  in  CNs  is  much  higher  than  in  superlattices.  To  explain  this  fact  let  us 
compare  the  mechanisms  of  the  nonlinear  conductivity  in  CNs  and  in  the  superlattices.  The 
quantum  superlattices  are  formed  by  alternating  plain  layers  of  different  semiconducting 
materials,  while  the  latteral  superlattices  consist  of  ID  chains  of  identical  and  identically 
coupled  GaAs/AlGaAs  quantum  dots  | ].  Both  these  structures  are  characterized  by  the 
dispersion  law  £(p)  — A|1  — cos(apz)|,  with  A as  the  overlapping  integral  and  a —2a/3. 
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Fig.  1. 


/— V curves  for  armchair  (left  panel)  and  metal  zigzag  (right  panel)  carbon  nanotubes. 


T = 287.5  K,  r = 3-  l(T12. 


Applying  the  method  described  above  to  this  dispersion  law,  one  can  obtain  the  Esaki-Tsu 
relation)  ] jz(Ez)  = crzz £z/(l  + £22t2)  instead  of  Eq.  (4).  Here  crzz  = \\mE.^o(djz/dEz) 
is  the  linear  conductivity.  Comparing  two  equations  for  the  current  density  we  can  con- 
clude that  the  specific  peculiarity  of  the  CNs  is  the  presence  of  the  high  Stark  components 
(summation  with  respect  to  r in  Eq.  (4)). 

It  has  been  shown  that  the  electron  motion  in  dc-field  can  be  described  as  the  os- 
cillations of  an  ensembles  of  effective  harmonic  oscillators  (Stark  components)  with  the 
frequencies  rQ.  The  full  current  is  a superposition  of  partial  currents  of  the  Stark  com- 
ponents. Their  electrical  field  strength  corresponding  to  maximum  current  of  r-th  Stark 
component  decreases  with  r as  r-1  while  our  calculations  shows  that  oscillator  strength 
of  these  component  decreases  slowly.  This  is  due  to  the  hexagonal  crystalline  structure  of 
CNs  reflected  in  dispersion  laws  (1)  and  (2).  The  number  the  unneglectable  components 
is  70-150  for  metal  and  200-300  for  semiconducting  CNs.  As  a result,  the  role  of  the 
high  Stark  components  in  CNs  is  essential  and  the  integral  nonlinearity  of  the  CNs  is  much 
higher  than  in  superlattices  [ ]. 

Note  that  the  NDC  provides  the  current  instability.  It  can  be  expected  that  simultane- 
ously applied  dc-  and  ac-fields  will  result  in  the  dynamic  electron  localization  (which  is  the 
nonlinear  phase  of  the  instability  ) and  in  the  2D  analogue  of  the  self-induced  transparency, 
like  it  takes  a place  in  the  semiconducting  superlattices.  The  above  mentioned  effects  are 
responsible  for  the  absolute  negative  conductivity  which  thus  is  predicted  to  be  exhibited 
in  CNs.  It  must  result  in  the  appearance  of  the  absolute  negative  conductivity  zones  and 
active  properties  of  CNs  providing  a potentiality  for  the  design  of  generative  nanodiodes 
in  microwave  and  infrared  ranges.  Such  a possibility  relates  to  both  single  CNs  and  CN 
ropes. 

The  NDC  effect  was  also  predicted  in  Ref.  [ ] for  the  CN-junctions  generated  in  the 
result  of  different  ways  of  doping  the  sides  of  the  tubes  or  via  the  applying  an  external 
potential  difference  between  them.  The  comparison  of  physical  mechanisms  of  the  NDC 
in  | ] and  [ ] is  actual.  In  | ] the  electric  field  is  localized  in  the  narrow  region  of  the 
junction  while  in  | j it  is  homogeniously  allocated  along  the  CN.  That  is  why  the  field 
strength  in  [ | is  much  higher  under  the  same  potential.  In  | ] the  NDC  was  obtained  for 
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the  junctions  based  on  metal  CNs  only,  while  [ ] also  contains  the  prediction  of  the  NDC 
in  semiconducting  CNs,  the  effect  being  much  weaker  for  them.  The  value  of  £^max)  for 
armchair  CNs  with  small  radius  is  approximately  half  as  great  than  that  for  metal  zigzag 
CNs  with  similar  radius.  This  result  was  obtained  in  both  papers.  Thus  one  can  come  to  a 
conclusion  that  physical  mechanisms  in  [ | and  [ ) are  close.  However,  it  is  necessary  to 
note  that  the  I—V  curves  in  [ ],  as  different  from  [ ],  are  unsymmetric  with  respect  to  the 
transformation  V — — V. 
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